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Functional Renormalisation Group (FRG) equations are constructed for a simple 
Yukawa-model with discrete chiral symmetry, including also the effect of a nonzero com¬ 
posite fermion background beyond the conventional scalar condensate. The evolution of 
the effective potential of the model, generically depending on two invariants, is explored 
with help of power series expansions. Systematic investigation of the effect of a class of 
irrelevant operators on the lower (stability) bound allows a non-perturbative extension 
of the maximal cut-off value consistent with any given mass of the scalar field. 

Keywords: functional renormalization group; Yukawa systems; stability bound; irrelevant 
operators 
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1. Introduction 

The experimentally found mass of the Higgs-field is light in the sense that it is 
very close to the lower edge of the perturbatively determined infrared window (for 
recent theoretical reviews, see EH). This situation resulted in an increased number 
of studies investigating the influence of quantum fluctuations on the stability of the 
Higgs effective potential. For the characterisation of this influence two conceptually 
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different interpretations were put forward. The first is motivated by the setup of 
Monte Carlo simulations of the ground state of field theories where one finds no 
signal for any instability in the continuum limit if the theory defined at the scale of 
the lattice constant (approximately the inverse momentum cutoff) is stable ffiEEEl 
Along this line of thought one arrives for any fixed value of the scalar mass at a 
maximal value of the cut-off beyond which no stable bare Lagrangian can exist with 
a fixed set of operator content. 

The alternative interpretation restricts the bare potential energy of the scalar 
field to the perturbatively renormalisable quartic power. One investigates the renor¬ 
malisation group equation (RGE) of the renormalized quartic coupling embedded 
into the standard model as a function of the renormalisation scale when the cut-off 
is sent to infinity | 8 | 9 | io | ii | i2 | i3 [ Qj^g encounters a negative value for this coupling 
upon fixing its ’’measured” value at some infrared scale, when the renormalisation 
scale goes beyond a certain, maximally allowed value. The order of magnitude of 
this ’’critical” value is the same as of the maximal cutoff value emerging from the 
first approach. 

The close agreement of the stability (lower) mass bounds obtained using the 
set of perturbatively renormalisable operators with various techniques ranging from 
perturbation theory, through semianalytic solution of the functional renormalisation 
group equations to the lattice simulations give considerable mutual support to the 
credibility of the results. 

Recently adding perturbatively irrelevant operators to the Higgs potential was 
shown to increase the maximal allowed cutoff considerabljlSIIEH It is not clear 
yet if there exist any absolute lower bound in this wider sense, since no systematic 
exploration of the coupling space has been conducted to date. 

In the present paper we extend the operator set involved in the determination 
of the lower mass bound with the ERG. Earlier approximate computations were 
restricted to a one-variable effective potential depending on the invariant p built 
exclusively from the scalar field. We shall consider the RG-evolution of the potential 
energy depending also on a second variable, the local invariant formed jointly with 
the scalar field and the fermions. The effect of this new invariant (called below the 
Yukawa-invariant, and denoted by I) will be investigated in the framework of the 
simplest scalar-fermion Yukawa theory 

Actually, we shall study the linear dependence on the Yukawa-invariant. Such 
an extension of the scalar-fermion action has been considered first in Ref llTI in the 
context of a quark-meson theory, where it emerges naturally from the RG-fiow. 
For strong enough scalar self-coupling similar dynamical generation might manifest 
itself in the present system. In first place, the point where the two-variable potential 
is expanded to linear order in the Yukawa-invariant will be varied. Apparently, this 
freedom leads to an almost negligible variation in the scalar mass bound at fixed 
cutoff. This feature justifies the omission of operators involving higher powers of 
I from the ansatz. The evolution of all couplings of the quantum action follows 
rather closely the rate determined by the perturbatively computed respective beta- 
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Scalar mass stability bound 3 

functions until only perturbatively relevant and marginal operators are considered. 
The lower bound characteristic for this case is reached when the initial value of the 
quartic scalar self-coupling is chosen zero. 

The situation changes spectacularly when one allows the coefficient in front of 
the Yukawa-invariant (the Yukawa-coupling) to depend on the scalar field invariant, 
introducing this way perturbatively irrelevant operators into the potential. In this 
case the evolution of the couplings determined by the RG-equations joins with a 
considerable delay the nearly perturbative regime after a non-perturbative ’’pre¬ 
history”. Initial data for the extra evolution are further optimised dynamically by 
running towards the ultraviolet from the optimal initial data found when restricting 
the initial quartic scalar coupling to zero. This running is stopped now when the 
model looses its stability through the generalised Yukawa coupling becoming neg¬ 
ative. This procedure is used to map the curve of the lowest mass estimate found 
for field independent Yukawa coupling onto a new curve. This way one extends the 
range of allowed cutoff momenta to values larger at least by one order of magnitude. 

The paper is organised as follows. In section 2 the model and the Ansatz for 
its scale dependent effective action are introduced. The set of RGB’s is constructed 
explicitly in presence of scalar and fermionic backgrounds in the formulation due 
to Wetterich [T8lT9](see also Ref i^Dll . The extra contribution reflecting the presence 
of the fermionic background is clearly identified. In section 3 the RGB’s are con¬ 
structed for the effective potential expanded in the Yukawa-invariant to linear order 
at an arbitrary point. Here we also describe the strategy for scalar mass estimates 
from the lowest non-trivial power-function ansatz for the potential. In section 4 the 
numerical solution will be presented in the Local Potential Approximation (LPA) of 
the Wetterich equation. We shall concentrate on the effects arising from the choice 
of the expansion point and from including the specific, perturbatively irrelevant 
piece into the Yukawa-coupling depending linearly on the scalar invariant. In sec¬ 
tion 5 the robustness of the effect of this operator is thoroughly investigated when 
one steps beyond LPA taking into account the anomalous dimensions of the fields 
(LPA') and also further higher powers are included into the field dependence of the 
Yukawa coupling and in the scalar potential. A summary and the conclusions of 
this investigation appear in Section 6. 


2. The model and the application of the Renormalisation Group 


The scalar-fermion Yukawa model in Euclidean space-time will be treated in the 
framework of the following effective action Ansatz defining the theory at scale k: 


Pfc = 




:^Z„k{dmcr)‘^ + Wk{p, I) 


p=-a'^, I = 


( 1 ) 

Here the last term represents the potential energy density of the model, which we 
shall expand to linear order in I: Wk « Uk{p) -I- hk(p)I- A new feature of the 
present investigation relative to Refill is the general p-dependence of the Yukawa- 
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type coupling. The model possesses a discrete 75 invariance: 

a{x) -a{x), tp^x) 75 V’(a;), ip{x) -tp{x)-f 5 . ( 2 ) 

T j r\ 

A non-zero homogeneous condensate at scale k given as Vk = <7k = y/2Z„kPk 

in the scalar field cr corresponds to the spontaneous symmetry breaking of the 
discrete 75 transformation. The condensate generates mass for the fermi-field: 

~ hk{pk)o'k, which itself varies with the scale. The physical condensate and 
the physical mass values are computed from the action T/j^oj at which the starting 
^k=A will land after running from the cutoff A to fc = 0: 

'i’ = zl()^ao, ml = 2hl{po)po, 

ml = U'M + 2poC/^'(po) + (3h'(po) + 2poh''{po))Io = U^^\po) + h^i^\p)h. (3) 


To ensure maximal analogy between our model and the top-Higgs sector of the 
Standard Model, in the discussion of the solution, below we fix the k = 0 values of 
the vacuum condensate to ?; = 246 GeV and the fermion mass to mp = 173 GeV. 

The RGEEIEHI as projected on § expresses the rate of variation of Tj. as a 
function of t = ln(fc/A) through an expression containing T^, the second functional 
derivative matrix of the effective action. It is represented as a 3x3 matrix, where 
the fields are represented by the transposed of the row-vector {a{x),tjj’^{x),'ip{x)) 
introducing the Gor’kov-Nambu representation of the fermi-field. One can transform 
this matrix into a block-diagonal form 1^, which leads for the second equality of 
RGE: 


5trfe[$] = -4striog(r(2) + R,) 

= -^4Trlog(r^"7 + ) + ^4Trlog(r^"j + ) 


-k-^tTrlog 


1 - 


.(2) 




-(2)_ 

<^bp(2) 


- Tf + 

'ipa 


n(2) 


R^ 


n(2) 

'ijj'R a 


U) 


Rk is the infrared regulator restricting the contribution to the functional trace Str 
onto the neighbourhood of the actual momentum scale k. The symbol dt refers to 
a t-derivative acting exclusively on Rk . Throughout this paper we shall employ the 
so-called optimised regulator The symbol denotes a hypermatrix in the 

doubled bispinor field basis where the condensed notation = {pjF and its 
transposed column-vector is used. Each term has a diagrammatic representation 
(see Figj^. The first two terms correspond to the free fermion and free boson 
contributions to the super-tracelog computed in the background Uk- The last term 
is the sum of the infinite series of 1-loop contributions containing an increasing 
number of alternating boson and fermion propagators. At each vertex of fermion-to- 
boson transformation an external fermion leg joins the loop, reflecting the composite 
fermionic background . 

In earlier papers the super-trace was evaluated on a background exclusively 
consisting of a nonzero u-condensate 1^1, where T^^^^ = = 0. In our recent 
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Fig. 1. Propagator sequences building the one-loop contributions to the tracelogs of . On the 
left a piece of a scalar loop (dashed curve), in the middle a fermion loop (continuous curve) appear 
both immersed into a scalar condensate, on the right a loop constructed of alternating sequence 
of fermion and scalar propagators immersed into a composite fermionic condensate is represented. 

publications clear arguments were put forward [13111 how the Grassmannian nature 
of the fermi fields is compatible with a coarse grained background where also nonzero 
n > 1 terms are present. Earlier this possibility was studied in a series of 
papers also by Aoki et al. p4 | 25 | 26 | 2 g ^ 


3. Renormalisation Group Equations in Local Potential 
Approximation (LPA) 

3.1. The potential energy 


For this calculation one chooses from the start space-time independent scalar and 
fermi backgrounds. On this background one has the following representations for 
the propagators using linear infrared regulator functions rB{k/p),rF{k/p)^^. 

= iZ^PF{p)i) - Pf = 1 + rF, TFijp) = O , 

( 5 ) 

and 

GaHp) = Z^Pb{p)p^ + mlj., Pb = 1 + tb, tbIp) = - 1^ 0 - 1^ . 

^ ^ ( 6 ) 

The masses are given with the expressions in (|^, but computed with the corre¬ 
sponding background fields and couplings defined at scale k. 

Below the notation will be used as a short-hand symbol for the momenta 
appearing in combinations ZaPB{q)<f or Z^Pp{q)q^, depending on whether one 
deals with the bosonic or fermionic propagator. One finds with a few algebraic steps 
the following expression for the rate of evolution of the potential energy density 
{h^k\pk) = hk{pk) + 2pkh'(pk)): 


^t[Uk{Pk) ^k{Pk)^k)] — 




log 




i<lR 


-km2)4 


log < 1 - 


hk{p){h^P {Pk)Y 2Ik 


<1r 


mf, 


Yn + ml 


( 7 ) 


In view of the general /fe-dependence of the right hand side, present also through 
nYr k-! right hand side is to be projected on the ansatz linear in figuring on 
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the left hand side. In fact, one can use an arbitrary function I^{pk) as expansion 
point: 

RHS(pfc,4) = RHS(pfc,4^(pfc)) + (4 - 4(p4)4-RHS(pfc,4(pfe)) + ..., (8) 

The projection has been realised in earlier investigations by an expansion around 
the origin 4 = O^Sl. If one would know the full two-variable potential W {p, I) then 
a natural choice would be 4 = Imin{Pk) determined by the condition dW/dl = 0. 
Such a procedure is analogous to the fermion-boson translation strategy of Reff251 
Keeping both terms displayed on the right hand side of Q and writing separately 
the RGE for the expression linear in 4 one has 


dtUkipk) = ^dt I <{-51og4^ 


‘^hkPk) 


+ log {q]i + ^ ^Ik{Pk)){q\ + ‘^h\pk) — 2hk{h\. ^)^4(Pfc)] } 






(4 


2hlpk) 


{ql + + h^k^Il{pk)){ql + 2hlpk) - 2hkihi^^rii 


■Ik, (9) 




^gf -h!i^\qR + 2hlpk) 

2 ‘ A {ql + + h^hl{pk)){ql + 2pkhl) - 2hk{h^pYIl ''' 


( 10 ) 

When one chooses p-independent Yukawa coupling and 4 = Oj ftio equations co¬ 
incide with the set of RGE’s of derived in Refill Still, the RGE of hk{p) alone is 
different depending on whether one defines it as F^J^^lSlor as dW{p, I)/dI as we 
do in the present study (see below). In both cases when one constructs a power 
series solution for hk(pk), then the most convenient is to expand it around pk,min, 
since the physical fermion mass is determined with ho(po,min)- 

For the solution of the RGE’s, they are rewritten in terms of dimensionless 
background field variables and couplings, rescaled by appropriate powers of the 
scale k. These quantities are distinguished by the lower index ’r’. The operation 
kdk is performed for the linear regulatoil^. The coupled equations of the potential 
and the Yukawa-coupling are as follows, after cancellation of 4 on both sides of the 
latter: 


dtUripr) + dUripr) + [2 - d - T]^)PrU!,{pr) 

5M^ {^ + + {1 + p'^)Mcr 

— — Z -9-1-77^;- 

[ 1 + il + pl){l + pl)-2K{hi^^)^lM\ 

+Vd -^- X 

[(l+/i2)(l+p2)_24(4l))2/.]2 

+ pl\) - 2hr{h'd^f‘){l+p^^)M^ - 24(hl^l)^(l -t p.^ 


42)4)M^1(,11) 
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dtK + {2- d - r]a)PrK + 2^4: - d - Per - 2prl,)hr 
-tu 

{hf\l + p%) - 2/i,(4i))2)(1 + 4)M, - 2hr{h^^^f{l + pl - 42)/«)m^](.12) 

Here the notations Vd = S'd/(d(27r)'^), V 4 ^ = (327r^)“^, = — 5* In Zo-, ??,/, = 

—9tlnZ^ and = 1 — Pa/{d + 2),M^ = 1 — p,p/{d + 1) are introduced, with 

Sd standing for the surface of the d-dimensional unit sphere. The boson (fermion) 
masses are rescaled by {Z^k'^) to give (m|)- 

In order the reader could appreciate the non-trivial numerical uniformity of the 
lo wer bounds presented below, the extra terms which would show up on the RHS 


of 


0 


when one defines h 


as are also given: 


+dt 


-dt 


2pU'") 


(g| + C/(l))2(g2 

1 (/i(l))3 


R w ‘2h‘^p) 




2/l2(/i(l))3 


C/(i))(g| + 2/i2p)2 






2h?p) 


(4 

q\ + C/(i) 


+ 


(<z| + c/W)2 


(13) 


3.2. Power series analysis for Ur and hr 

The simplest search for the solutions makes use of power series of Ur and hr and 
sets = Za- = 1. One finds nontrivial results already with a quadratic expres¬ 
sion for Ur which corresponds to keeping in the effective action only perturbatively 
relevant-l-marginal operators. The flexibility of the lower bound under some mod¬ 
ulation of the Yukawa coupling will be investigated with an ansatz linear in pr'. 
hr = hf) + hiPr- In the symmetric phase one has 


Ur{pr') — ^IkPr 2^‘^^Pry 


Pp — 2(/l0 + hiPrYPr, P^r — Aife -|- iX2kPr + 3/li/® 


(14) 


(Aifc > 0), while the solution in the broken symmetry phase is parametrized as 

Ur{^P — Pr Pr.min) ~ 2 ’ Pa ~ ‘^X2kPr,min ~h 3X2k^P ~h 3h\I^, 

Pp — Pr,min + 2(h^ + 2hhiPr^min) ^P + ‘^hiiJliPr^min + 2h)Ap^ + 2h'\Ap\fb) 

with Pr,min, the rescaled minimum of Ur, {Ur{pr,min) = 0) and h = hg + hipr^min- 
The above mentioned two cases for the choice of Ir{pr) were compared. After 
expanding the right hand side of RGE in power series of pr or Ap and equating 
the corresponding terms one arrives at a set of equations for the couplings A„fc 
and/or the expectation value Pr,min- The consistency of the two sides are ensured 


by keeping terms on the RHS of (111 to quadratic, while on the RHS of (12) to first 
order in pr or Ap. 

The discussion of the influence of the truncation at higher (perturbatively irrel¬ 
evant) couplings of Ur on the stability bound requires some care since truncations 
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at odd powers lead to seemingly unstable potentials with the coefficients in front of 
the highest power converging to zero from the negative side when k —>■ oEUH Simple 
argument shows that this is an artefact of the power series truncation (for a recent 
analysis in 0(N) models see R ef IH(7|l . Actually, resummation of any power series 
expansion (for instance with help of Fade approximants should reconstruct the 

d 

correct asymptotic behaviour of Ur a t the fixed point, which is ^ pr ~^, as it can be 
seen from asymptotic analysis of ( |ll| . Therefore Ur = f{pr)pr ~^, where f{pr) —t 1 
as pr —t oo. As a consequence, the coefficients of the power series expansion of f{pr) 
have alternating sign. This artefact appears explicitly when in section 5.3 the term 
XskP^/S is included into Ur- The infrared values of the physical quantities are not 
affected by the seemingly unstable behaviour of the odd truncations. 

The effect of a different class of irrelevant operators on the stability bound can 
be studied assuming a non-trivial p-dependence of the Yukawa coupling. In the 
above explicitly described simplest non-trivial example of involving an irrelevant 
operator one might test its effect on the stability bound by starting the RG-flow 
with /io(fc = A) « 0 and /ii(fc = A) > 0 and tuning hQ{k = 0) + hi{k = 0)po,min to 
the value determined by the fermion mass. 

4. Main result: non-perturbative extension of the stability range 

through an irrelevant Yukawa-type coupling 

First, we summarize the findings of ReffHl concerning the edge of the stability 
range without perturbatively irrelevant operators. In agreement with the one-loop 
perturbative analysis the lower bound from the RGB’s of the system stems from a 
set of initial data characterised by A 2 (A) = 0. The quartic coupling evolves from 
here essentially logarithmically (see the dashed line of Figj^ until it is saturated 
with its infrared (k = 0) value at around —t > 10, which determines the physical 
scalar mass. The resulting curves of the cut-off dependence of the lower scalar mass 
bound computed with one of the two values of differ from the results of Ref[T4] 
only slightly as one can see looking at the two nearby lying curves corresponding 
to h(p) = ho in the left hand side of Figj^ The difference is the consequence of the 
different definition of hk, mentioned above and leading to different RGB’s beyond 
the leading one-loop beta-function. The choice of the expansion point has even 
less (fully negligible) effect: the respective curves coincide within .1%, therefore 
cannot be resolved in the figure. 

Concerning the influence of including perturbatively irrelevant operators we do 
not discuss in this paper standalone the effect of switching on higher powers of pr 
in Ur{pr) which was investigated earlier in Refs fHl [TSl We concentrate on allow¬ 
ing linear p-dependence in the Yukawa-coupling {hi ^ 0) which by itself leads to 
interesting and potentially important consequences. 

The non-trivial running of hi has the effect of diminishing the starting slope of 
the A 2 -evolution, when one starts from A 2 (A) = 0 towards the infrared (negative 
t-values). The starting slope gradually decreases with increasing hi{A) (see the con- 
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Fig. 2. The evolution of A 2 towards lower scales {t < 0) for A = 10^ GeV starting with \ 2 {k = 
A) = 0. The dashed line corresponds to hi = 0. The continuous line represents the evolution for 
hi (A) = hi^rnax when A 2 starts with zero slope. In this case one can smoothly evolve the couplings 
into the region t > 0. 


lf=0 and /7*=rg)p^(0,0,0): 

h^=Q 

- 

h-\*0 + UV extension 

if=if{Eow\y. 

- /}1=0 

- /)1*0 

- + UV extension 

104 105 106 107 108 109 

A 

Gev 

Fig. 3. The upper left lying two continuous curves represent the lower bounds mo-,min (A) obtained 
for hi = 0 and for hi = hi^max, starting the RG-evolution from the scale k = A and using the 
expansion point I^{EoM). The considerably lower bound curve in the middle of the figure is 
obtained by extending the allowed momentum range dynamically to A' = Aexp(£macc)- Also the 
three corresponding bounds arising with the RGE of h{p) derived by Gies et al. and using 7^ = 0 
are presented (dashed curves). 



tinuous line in Figj^. It rectifies to its original value with a slight delay, reducing 
the value of the lower bound belonging to the chosen cutoff by a small amount. The 
largest acceptable hi (A) value (denoted as hi „iax) corresponds to A 2 starting hori¬ 
zontally, since beyond this value A 2 (fc) would (temporarily) cross into the unstable 
regime. The corresponding slightly lower bounds lie below, but close to the bounds 
obtained with hr{pr) = h-o- 

For an essential decrease in the lower bound the following observation is of cen¬ 
tral importance. In the zero slope case one can continuously integrate the RGB’s 
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t 


Fig. 4. The zero of ho arising when one evolves it into the t > 0 region stops the dynamical 
extension of the cutoff theory 


also upwards {t > 0). One finds a steeply increasing curve for X 2 {t > 0) also dis¬ 
played in Figj^ This segment of the evolution of A 2 is controlled by the competition 
of two terms on the right hand side of its RGE. The term ~ —Hq comes from the 
fermion quadrangle and tends to diminish A 2 . The other term ~ +hQhi increases 
A 2 and arises from projecting the contribution of the fermion bubble on the ^ 
term of Ur- It turns out that first stays constant, then for t>l steeply decreases 
(see Fig|4|. On the other hand hi displays ~ scaling behavior, therefore dom¬ 
inates for t > 0. As a consequence A 2 is expected to scale ~ which in fact is 
observed. Such a change of the scaling regime of A 2 when turning from t < 0 to 
t > 0 is not exceptional. A qualitatively similar change also happens around the 
point A 2 = /3a 2 = 0i ’’forced” on the evolution of A 2 by coupling the Standard Model 
to an asymptotically safe quantum gravity theorjJ^. 

The evolution for t > 0 eventually results in a strongly coupled theory with 
nonzero scalar condensate. The value of the condensate is 0{A' = Aexp(t)) and with 
it both the scalar and the fermion masses are of the order of A'(t): the scale hierarchy 
disappears towards the ultraviolet. The integration is stopped at the effective cutoff 
when the steeply decreasing value of ho (see Fig|^ reaches zero (a negative 
Yukawa-coupling (for any value of p) would destabilize the theory). In our experience 
no Landau-type (triviality) singularity showed up in A 2 until this point was reached. 

With the above procedure one determines dynamically an effective cut-off value 
A'max with the corresponding couplings (among them ho{A'^^^) = 0) where one 
actually can start the RG-evolution. One arrives at the same infrared data as if 
one would start at t = 0, but now the low-scale (physical) data are hierarchically 
separated from the initial data set. The new piece of the RG-trajectory, which one 
can call as strongly coupled ’’prehistory”, causes a substantial delay in the start of 
the logarithmic variation of A 2 , and results in a significant extension of the physical 
momentum range (see the right hand curve in Figj^. The RGE-variant of ReffHl 
can also be treated this way and the emerging new lower mass bound (which appears 
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first in the present paper) is also displayed in Figj^ 

Before closing this section we extend our analysis to the A 2 (A) > 0 case. In 
the light of the above discussion it is obvious that there exists also in this case 
a specific huuned choice for which = 0. More than this is true: a whole 

interval himin < huuned < himax can be stretched where the ’’partial” fixed point 
and the change in the scaling of A 2 occurs, just shifted to slightly higher or lower 
scales (kpp). The corresponding value of \ 2 {kFp) is below A 2 (A), but is required 
to stay non-negative. This observation makes it clear, that a possibility for smooth 
UV-extension works for the full allowed infrared range of the Higgs-masses without 
any need to fine-tune hi. The need for some sort of fine-tuning appears only when 
one aims at finding the lower mass-bound. 

5. Robustness of the UV-extension within a wider operator set 

In the previous section it was demonstrated that the effect of hi on the evolution of 
A 2 leads to the emergence of a ’’partial” fixed point where fi\^ = 0. Starting from 
this point a smooth continuous RG-extension of the validity range of the theory is 
possible towards the ultraviolet. In this section we present a number of evidences 
for the robust nature of the appearence of the ’’partial” fixed point of A 2 also upon 
taking into account the anomalous dimensions of the fields (LPA’) and including 
further perturbatively irrelevant operators into the RGE analysis. 


5.1. The impact of anomalous dimensions 


The two RGE’s characterising the evolution of the effective potential are completed 
by the algebraic equations which determine The projection of Q on the 

coefficients of the kinetic terms by consistency requires the evaluation of the right 
hand side at some specific value of the background fields. Two expansion points are 
tested in this paper. Either one chooses /^ = 0 or we use the point solving the field 
equation at scale k: 

f^k\i’'^)k + crkUk{pk) = 0 ^ h^^'> Il{EoM) + 2pk%{pk) = 0, (16) 

both imply = 0. This choice leads back to the RG-equations for and 

determined by 1^3 


Ptp — 2.V(ih^ 


(1-f/i2)2(l + ^2) y 


1 - 


her 


(17) 


= 2Vd 


{iU'f + 2prU'f'f 
2hl _ 2 


m 




1 + m; 


I 

1 


1 1 - w 

2 d- 2 


■ ( 18 ) 
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t 


Fig. 5. The RG-evolution of the anomalous dimensions rja- and 77 ^ (A 2 (A) = 0,A = lO^GeV. 
Initial values of other couplings are tuned as indicated in the main text.) 


Here is the scaled scalar mass with no correction from I^{pr)l Since the RG- 
evolution of the model ends in the broken symmetry phase where both fields are 
massive and therefore ~ these equations imply that —>■ 0 when 

k ^ 0. The anomalous dimensions actually might influence the evolution of Ur and 
hr only at intermediate scales. The ry-evolution is exemplified in Figj^where for A = 
lO^GeV and A 2 (A) = 0 the scaled (dimensionless) values of the other three couplings 
are tuned to the values Ai(A) = 0.007188, /io(A) = 0.7546, hi{A) = 0.4295, in order 
to have the physical values of and v in addition to arranging (^) = 0 . 

From and ( |18[ ) one can express the anomalous dimensions algebraically and 
substituting them into Eqs. (11 1 and (12) one proceeds to the solution of the result¬ 
ing equations with the same algorithm as applied when = Zp = 1 (this is the 
LPA' variant of the local potential approximation (LPA)). Our general experience 
is that the anomalous dimensions stay very small (see Fig[5]) , therefore the lower 
mass-bounds obtained in LPA remain unchanged even quantitatively. 


5.2. The impact of the strength of the operator 

This operator takes over the stability control of the potential Ur along the p- 
direction from the ~ p^ term and lowers considerably the lower Higgs-mass bound 
as was demonstrated in Pefs lMlITHl Here we consider a similar situtation where the 
RG-evolution towards the infrared starts with A 2 (A) < 0 ,A 3 (A) > 0,A = lO^GeV. 
As one sees in Figj^the freedom to tune Ai(A) = 0.00736, ho(A) = 0.756,/ii(A) = 
1.1 allows, in addition to fixing the physical values of v and mp, also in this case 
to tune = A) = 0. This in turn opens the way to the UV extension of the 

theory up to some A' > A. A 3 has a diminishing effect on the rate of A 2 and eventu¬ 
ally compensates at some positive t-value the increasing tendency experienced for 
A 3 = 0. This compensation is followed by a fast crossing of A 2 to negative values as 
one sees in Fig|^ However, in this case the negative value of A 2 does not prompt for 
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A2(A3 = 0) 

MM > 0 ) 
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holM > 0) 
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20 


Fig. 6 . Illustration of the effect of switching on A 3 on the RG-evolution of A 2 and ho (A 2 (A) = 
—0.05, A 3 (A) = 3, a = lO^GeV. Initial values of other couplings are indicated in the main text). 

any stability consideration, since A 3 > 0. The UV-evolution is stopped by the zero 
of Hq which happens slightly earlier than for A 3 = 0. One also notes the significant 
decrease in the infrared value of X 2 {k = 0 ) relative to the A 3 (A) = 0 , A 2 (A) = 0 
case. Therefore the effect pointed out in ReffHl combines here with the effect of 
p-dependent Yukawa-coupling, to lower further the lower Higgs-bound obtained by 
applying the two irrelevant operators separately. 



i 


5.3. The impact of the strength of the operator p^I 


The RGE for the coefficient /12 can be deduced by projecting Eq.(12) on this oper¬ 


ator. One tunes also here the couplings Ai(A) = 0.007168, ho(A) = 0.756, hi (A) = 
0.4153 to ensure m^,v, = 0. 

It is important to recognize that h 2 does not appear directly in the RGE of A 2 
(one easily figures out that from the fermion-bubble, the lowest power of p where 
h 2 contributes to the RGE is p^). On the other hand h 2 shows up in front of a 
tadpole type contribution in the RGE of hi. Towards the infrared for moderate 
( 0 ( 1 )) values of ^2 there is no observable change in A 2 (fc = 0 ). 

The more important question refers to the effect of h 2 on the UV (t > 0) 
extension of hi. There the increase of hi is slowed down by its effect, and eventually 
after reaching a maximum hi starts to decrease. Since hg behaves the same way as 
for h 2 = 0, one reaches a t-value where the term ^ —/iq compensates ^ +hohi. In 
this point A 2 reaches its maximum followed by a fast drop. The UV-extension is 
stopped by the zero of A 2 or that of hg, occuring first. This evolution of A 2 and of 
hi is illustrated in Figj^by comparing the RG-evolutions with ^2 = 0 and /i 2 ^ 0. 
Increasing / 12 (A) pushes the zero of A 2 closer to the location of the other zero, 
the origin (for the present set of parameters the inflection happens for /i 2 ~ 65). 
Anyhow, there is a natural (and rather wide) range of / 12 (A) ~ 0{1), where the 
shift of the curve of the lower mass-bound is the same as for /12 = 0 . 
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Mhi = 0 ) 


Mh2 > 0 ) 
> 2 > 0 ) 


Fig. 7. Illustration of the effect of switching on h 2 on the RG-evolution of A 2 and hi (A 2 (A) = 
0,h2(A) = 0.5, A = lO^GeV. Initial values of other couplings are tuned as indicated in the main 
text.) 


Based upon these three specific cases, we conjecture the generic nature of the 
effect of the operator hi pi on the evolution of A 2 resulting in a ’’partial” fixed point, 
where = 0. Starting from this point one is allowed to continue the RG-evolution 
of A 2 towards the ultraviolet. The Landau-singularity of A 2 is apparently avoided, 
the range of extension is limited by the instability caused by other operators. 

6. Conclusions 

The motivation for the presented investigation was to explore the RG-evolution of 
the two-variable effective potential W {p, I) depending also on the fermionic back¬ 
ground through I. We have derived RGE’s for its expansion in I to linear order, 
and investigated the effect of the choice of the expansion point /*. 

In this framework it is natural to investigate the case of a genuine p-dependent 
coefficient in front of / in the quantum action. The emergence of the perturbatively 
irrelevant operators is quite natural in the FRG-paradigm. In the absence of a self- 
consistent UV-completion (governed by a new fixed point), the extended effective 
action simply describes a cut-off theory. 

We have performed a search for the lower scalar mass bound by tuning the 
dimension-6 coupling hi, and optimised the extension of the stability range of the 
theory along this direction. It was demonstrated that the UV-extension of the range 
of validity of the field theory is made possible through the existence of a ’’partial” 
fixed point {j3x^ {kpp) = 0) for A 2 (A) > 0 showing up for some scale kpp by choosing 
the value of hi from a finite interval. These findings were independently confirmed 
by the Jena-groufP^. In an extended investigation we also have checked that the 
location of this interval is slightly shifted when other perturbatively irrelevant op¬ 
erators are also switched on or the anomalous dimensions of the scalar and fermion 
fields are taken into account. 
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The theory arrived at by the UV-extension, at high energy is a strongly coupled 
theory with a ground state of broken discrete chiral symmetry. All masses are of the 
order of the cut-off, there is no mass hirarchy in the ultraviolet. The RG-evolution 
first leads to vanishing /3-function of the quartic scalar coupling A 2 at an interme¬ 
diate scale, wherefrom the theory closely follows the perturbative evolution to the 
final infrared mass values hierarchically separated from the cutoff. This construction 
is hopefully applicable also in other models involving Yukawa-type interaction. 

All investigations of the quantum stability of the present model performed to 
this date point to the possibility of a non-perturbative extension of the allowed mo¬ 
mentum range with the application of perturbatively irrelevant operators. One way 
to further check the persistence of this feature is to continue the gradual expansion 
of the coupling space along the lines presented in section 5. A more comprehensive 
approach is to attempt the complete determination of the scale dependent two- 
variable potential W{p, I) with the correct large p and I asymptotics as has been 
realised very recently for the one-variable U (p) functioiP^. 
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